ARPES Spectra of the Hubbard model 
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We discuss spectra calculated for the 2D Hubbard model in the intermediate coupling regime with 
the dynamical cluster approximation, which is a non-perturbative approach. We find a crossover 
from a normal Fermi liquid with a Fermi surface closed around the Brillouin zone center at large 
doping to a non- Fermi liquid for small doping. The crossover is signalled by a splitting of the Fermi 
surface around the X point of the 2D Brillouin zone, which eventually leads to a hole-like Fermi 
surface closed around the point M = (7r,7r). The topology of the Fermi surface at low doping 
indicates a violation of Luttinger's theorem. We discuss different ways of presenting the spectral 
data to extract information about the Fermi surface. A comparison to recent experiments will be 
presented. 



Introduction 

The rich phenomenology of high- Tc superconductorsEI 
has stimulated strong experimental and theoretical inter- 
est in the field of strongly correlated electron systems. 
Apart from the anomalously high transition tempera- 
tures, these compounds are also of interest due to their 
unusual normal state properties. Most of these anoma- 
lous properties are found in spectra and transport quan- 
tities, i.e. are intimately linked to the dynamics of the 
electronic degrees of freedom. Thus, much of the exper- 
imental and theoretical effort has concentrated on the 
development of an understanding of the single-particle 
dynamics. Among the fundamental and controversial 
questions are whether the cuprates can be described as a 
Fermi liquid or not and what shape and volume a possible 
Fermi surface will have. 

In this connection, one of the most informative exper- 
imental probes of the cuprates has become the Angle- 
Resolved Photoemission Spectroscopy (ARPES). The de- 
velopment in this field, again largely stimulated by the 
interest in the physics of high-Tc, has led to aXtemendous 
increase of the angular and energy resolutionElcl. ARPES 
is now able to access the low-energy behavior of single- 
particle spectra and especially the shape and topology of 
the cuprate Fermi surface.- This has lead to the discov- 
ery of the shadow bandEElu andnthe pseudo-gap forma- 
tion in the underdoped cupratesEJ. Recently, the single- 
particle self energy was extracted from the ARPES data, 
showing extremely interesting behavior especially close 
to (7r/2,7r/2) on the Fermi surfaced. 

A large number of the more recent ARPES exper- 
iments have concentrated on Fermi surface mapping. 
Some of these experiments seem to indicate that, at least 
for LSCO, the Fermi surface switches from being hole- 
like, centered at M = (7r,7r), for low doping to being 
electron-like, centered at the zone center F = (0, 0) for 
high|-|drtping with a volume consistent with the doping 
leveB'QQ. In other results, Pjacticularly for bilayer com- 
pounds like YBCO or Bi2212BEl, the Fermi surface seems 
to remain hole-like independent of the doping level. Es- 



pecially for those latter compounds, further complica- 
tions in the experiments arise from superstructures due 
to umklapp scattering and \hp. possibility of bilayer split- 
tings in YBCO and Bi2212l3tl. Furthermore, in the un- 
derdoped regime one can observe additional changes in 
the Fermi surface topology that are frequently brought 
in connection with the possibility of stripe formation in 

Lscoa 

The presence of a large Fermi surface has been taken 
as a validation of Luttinger's theorem; however, other re- 
sults find that the Fermi surface, volume, at least in the 
underdoped regime, is too smalH, which would point to- 
wards a violation of Luttinger's theorem. Furthermore, 
recent experiments seem to indicate that the low doping 
Fermi surface near X = (tt, 0) actually bifrpxates into two 
parts, one electron-like and one hole-likeM. This split- 
ting has-been interpreted in terms of strong interlayer 
couplingEI, but could equally well result from shadow 
Fermi surface formation due to coupling of the electrons 
to strppg antiferromagnetic fluctuations around the X 
pointsu. 

Thus, a consistent experimental picture concerning 
both shape and volume of the Fermi surface is at present 
not available and a theoretical investigation of the generic 
features to be expected based on a model calculation is 
necessary. 

Early in the theoretical investigation of the high-Tc 
cuprates it was realized that the 2D Hubbard model 
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in the intermediate coupling regime, or closely related 
modelsrlike the t-J model probably capture the essential 
physicalil. In the wake of this conjecture, a huge effort 
has been directed to the study of these modelsE3. There 
is now a general consensus that the appropriate param- 
eter regime for the cuprates is the intermediate coupling 
regime where the Coulomb parameter U is roughly equal 
to the bandwidth. However, this is the most compli- 
cated regime of the model since both weak and strong 
coupling perturbative approaches fail. Exact diagonal- 
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ization of small clustersta suffers from strong finite-size 
effects, often ruling out the reliable extraction of low- 
energy dynamics. Conventional Quantum Monte Carlo 
for finite sized systems suffers from a severe minus sign 
problem in this parameter regime. The resulting data is 
of insufficient quality to allow for reliable calculations of 
dynamic quantities at low enough temperatures. High- 
temperature series has provided some of the most infor- 
mative results for the Fermi surface topology, but it does 
not yield spectra, aiid, so far, only results for the t-J 
model are availableEj. 

Thus, also from a theoretical point of view, the ques- 
tion of whether a Fermi surface does actually exist and 
what its topology is still is a matter of debate. In that 
connection it is of special interest that some experiments 
indicate a violation ofpLuttinger's theorem; if true, any 
theory, such as FLEXEJ, based on a weak-coupling ex- 
pansion around the non-interacting limit would be inad- 
equate. 

Thus a treatment within a non-perturbative scheme 
clearly is desirable. In this paper we therefore use 
the recently dp'^celpped dynamical cluster approximation 
(DCA)yMO'BE3 to study the low-energy behavior of 
the 2D Hubbard model in the intermediate coupling 
regime with nearest-neighbor hopping t and on-site cor- 
relation U equal to the band width W. The DCA sys- 
tematically incorporates non-local corrections to local ap- 
proximations like the dynamical mean field, by mapping 
the lattice onto a self-consistently embedded cluster. We 
solve the cluster problem using a combination of quantum 
Monte Carlo (QMC) and the maximum entropy method 
to obtain dynamics. This technique produces results 
in the tkermodynamic limit and has a mild minus-sign 
problemtj. 

The paper is organized as follows. The next section 
contains a brief introduction to the DCA. The numerical 
results will be presented in the third section followed by 
a discussion and summary. 



can be varied systematically from short ranged to long 
ranged by increasing the number of coarse-graining cells. 
The DCA collapses to the DMFT if one represents the 
Brillouin zone by one cell only, thus setting the charac- 
teristic length scale to zero. 

By construction, the DCA preserves the translational 
and point group symmetry of the lattice. Comparisons to 
otheij-ej&tensions of the DMFT developed during the past 



yearSifj either show that these are identical to theX)CA 
or converge more slowly as function of cluster sizeE3. 

For the impurity problem of the DMFT a large set 
of reliable numerical. techniques has been developed over 
the past ten yearscJCJEj. We have employed Quantum 
Monte Carlo-ZQMC) and the Non-Crossing Approxima- 
tion (NCA)ll3 as non-perturbative (in U) methods to 
solve the cluster problem of the DCA. For cluster sizes 
larger than Nc — 4: however, only the QMC technique is 
presently available for the DCA. The cluster pmblem is 
solved using the Hirsch-Fye impurity palgorithmO modi- 
fied to simulate an embedded clusteiO. Note that this 
problem is computationally much more difficult than that 
encountered in either the DMFA or in finite-sized simu- 
lations, since the block diagonal structure in space-time 
occurring in conventional finite-system simulations is not 
present here. This increase in computation time is, how- 
ever, partially compensated by a rather mild minus-sign 
problem, even for camparatively large values of U and 
small temperaturesE3. From the QMC data, the spectra 
are obtained by^aiialytic continuation with the maximum 
entropy methocuJ. Finally, the self energy is interpolated 
on to the full Brillouin zone using Akima splines, which 
is a sensible step as long as the assumption of a slow vari- 
ation in k-space is valid. Note that it is very important 
to interpolate irreducible quantities like the self energy 
and not for example the cluster Green function itself. 



Results 



Formalism 

A detailed discussion of the DCA jfoMMjalisffla-jwas al- 
ready given in previous publicational3-t3liZl'll3lij. The 
main assumption underlying the DCA is that the single- 
particle self-energy z) is a slowly varying function of 
the momentum k and can be approximated by a constant 
within each of a set of cells centered at a correspond- 
ing set of momenta K in the first Brillouin zonec£l. The 
single-particle lattice Green functions are then coarse- 
grained or averaged within these cells, and used to cal- 
culate the lattice self energy and other irreducible quan- 
tities. Within this approximation, one can set up a self- 
consistency cycle similar tp-the one in the dynamical 
mean- field theory (DMFT)E3. However, in contrast to 
the DMFT, where only local correlations are taken into 
account, the DCA includes non-local dynamical correla- 
tions. The length scales of these non-local correlations 



For a proper description of the Cu02 planes of the 
high- Tc cuprates within the Hubbard model (|l|) it is gen- 
erally accepted that the tight-binding dispersion has the 
form 

ig= -2t (cosika;) + cos{ky)) (2) 

— 4t' cos{kx) cos{ky) 

with a nearest neighbor hopping amplitude t > and 
a next-nearest neighbor hopping amplitude t', which in 
principle can have any sign. From bandstructure calcu- 
lations and the general form of the measured Fermi sur- 
face, especially in the pvcrdoped regime, conventionally 
a negative t' is inferredn3. Such a negative t' would natu- 
rally lead to a Fermi surface closed around the Brillouin 
zone corner M. The interesting question, however, is 
whether the Fermi surface closed around the M point ob- 
served experimentally in the low-doping region is a mere 
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FIG. 1: (a)-(c) The single-particle spectrum A(k,t<;) for 
U = 2eV, T = ^eV, & = 0.05, t± = and iVc = 16 along 
certain high symmetry directions (spectra for different k are 
shifted along the y-axis). The colored lines in Figs, (a) and 
(c) indicate the spectra which cross the Fermi energy with a 
peak closest to lj = 0. In (b), no such peak is found which 
crosses the Fermi energy, (d) the imaginary part of the self 
energy versus frequency at the Fermi surface crossing found 
in (a) and (c). 



bandstructure effect or induced by correlations. Weak 
coupling treatments of the 2D Hubbard model indicate 
that such a change of the shape, but not the volume, 
of the Fefrn i surface due to the interactions is indeed 
possibleEJo. Thus, to obtain insight into the effects of 
correlations on the structure of the Fermi surface and dis- 
tinguish them from pure bandstructure effects, we con- 
centrate on the case i' = in this paper. 

In the following we set t = l/4eV in accordance with 
typical values extracted from the experiments and band- 
structure and choose U — W = 2eV. This value of U is 
sufficiently large that icf^Nc > 4 a Mott gap is present 
in the half- filled modelEJ. We performed our simula- 
tions at a range of temperatures, but will present results 
for T — 0.033eV only, which is roughly room tempera- 
ture. A pseudogap due to short-ranged spin correlations 
is also present in the weakly doped model for slightly 
lower temperaturestj. 

The single-particle spectra for certain high symmetry 
directions are plotted in Figs. H and || for n = 0.95 and 
n = 0.80, respectively. We use the standard convention 
to identify the high symmetry points in the zone, F = 
(0,0), M = (7r,7r) and X = (7r,0). For n = 0.95 the 
peak in the spectrum crosses the Fermi energy along the 
F M and M X directions, while for n = 0.80 the 
second crossing appears along X ^ T. The imaginary 



FIG. 2: (a)-(c) The single-particle spectrum A(k,di;) S = 
0.20 along certain high symmetry directions. Other parame- 
ters as in Fig. ^. The colored lines in Figs, (a) and (b) indicate 
the spectra which cross the Fermi energy with a peak closest 
to = 0. In (c), no such peak is found which crosses the 
Fermi energy, (d) the imaginary part of the self energy versus 
frequency at the Fermi surface crossing found in (a) and (b). 



part of the self energy at these crossing points is plotted 
versus frequency in panel (d) of Figs. |] and |^. 

One very interesting feature of the spectra at low dop- 
ing. Fig. |l|, is that the peak near (7r/2,7r/2) broadens 
dramatically before crossing the Fermi energy. Near X, 
on the other hand, one does not observe any dramatic 
change in the spectrum when crossing the Fermi en- 
ergy. This indicates that near (7r/2,7r/2) hole-like quasi- 
particle (QP) excitations with k < kp appear to have 
longer lifetimes than electronic excitations with k > kp. 
This asymmetry between particles and holes near the 
Fermi surface is a strong indication of non-Fermi liq- 
uid (NFL) behavior, at least along the T-M direction. 
Note that this is at least qualitatively in agreement with 
experimental ARPES spectraB'Ell. Although these exper- 
iments on LSCO hjad to be performed below Tc in the 
underdoped sampleE^, the authors infer from the general 
behavior of the spectra, that a peak crossing the Fermi 
energy can be found along AI X, while a rather broad 
structure is seen along F ^ M around (7r/2,7r/2). It 
would clearly be interesting to have experimental ARPES 
data above Tc available for a more founded comparison. 

It is also quite instructive to look at the imaginary part 
of the self-energy at the k-points where the peak in the 
spectrum crosses the Fermi energy (Fig. |l|(d)). In partic- 
ular at the crossing point kf, close to (tt, 0), 3ml](kb,a;) 
shows a striking asymmetry in the low-frequency regime 
as compared to 3mI](kQ,w), where ka is the crossing 



point close to (tt, 0). In fact, SJmI](kf,,w) starts to de- 
velop an additional feature which eventually leads to 
the formation of ajpseudo gap in the spectra along the 
M — > X directionEj. In addition, one observes a rather 
large residual scattering rate for both momenta ka and 
kb. This can either be taken as further evidence for 
NFL behavior or as signal-ior the occurrence of a new 
very small low-energy scaleE3. Obviously, we cannot de- 
cide this question on the basis of the present data, but 
would have to look at much lower temperatures. Unfor- 
tunately, this is not possible at present. Note that in 
the low-energy regime the self energy displays significant 
fc-dependence. This clearly renders theories based on a 
local approximation like the DMFA inadequate at least 
for small doping. 

At higher doping. Fig. ||, the peaks in the spectrum 
close to the Fermi energy are far sharper. Here it makes 
sense to speak of a conventional Fermi liquid and quasi 
particles again. As already mentioned, the Fermi en- 
ergy crossings can be found along F M and X ^ T. 
Again, this is in qualitative accordance spih ARPES 
experiments for strongly overdoped LSCOBEiJ, although 
these experiments still find rather broad structures even 
in heavily overdoped samples. In addition, there is no 
evidence for particle-hole asymmetry in our data. Es- 
pecially at (7r/2,7r/2) the structure crossing the Fermi 
level appears to be rather symmetric with respect to the 
crossing point. 

The imaginary part of the self energy, shown in 
Fig. ^(d), has a broad maximum at a; = with a very 
small residual scattering rate and changes little as k 
moves along the Fermi surface. This weak dependence 
on k is an indication that approximations like the DMFA 
should be accurate here, i.e. that there is little effect 
of non-local correlations. All indications are that for 
this doping regime standard Fermi-liquid behavior has 
returned. 

More evidence for NFL behavior can be seen in the 
shape of the Fermi surface. Theoretically, the Fermi 
surface can be defined in two different ways. First, the 
gradient of the electronic distribution function, |Vn(k)|, 
has a maximum at the Fermi surface. From a compu- 
tational point of view this quantity is very convenient, 
because it does not require the calculation of dynamical 
properties. However, especially for a comparison with 
experimental Fermi surface mappings based on ARPES 
experiments, the approach via |V7i(k)| probably is not 
adequate, mainly due to the unknown influence of ma- 
trix elements in the experimental spectraEI. 

An analysis of the Fermi surface for the t-J model 
based on a study of |Vn(k)| has beenj^erformed recently 
within a high-temperature expansionES and been consid- 
ered as clear evidence for a violation of Luttinger's theo- 
rem and the possible formation of a non-Fermi liquid at 
small doping. Our results for that quantity are collected 
in Figs. ^ and ^, where | Vn(k)| for the upper right quad- 
rant of the first Brillouin zone is shown in a density plot 
for n = 0.95 and n = 0.8, respectively. Regions of large 
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FIG. 3: Gradient of the distribution function |Vn(k)| for U — 
2eV, T = ^eV, S = 0.05, = and = 16. Shown is 
the upper right quadrant of the first BZ. The violet (red), 
represents regions of high (low) values. The blue and violet 
regions map out the Fermi surface. The non-interacting Fermi 
surface is shown by the black line. The interacting Fermi 
surface is rather hard to define. Note the strong bifurcation 
around the X points. 



values of |Vn(k)| are colored in blue and violet, regions 
of small values in red. For comparison the Fermi surface 
for the non- interacting system is included (black line). 
For small doping. Fig. ^, |Vn(k)| gives a rather broad 
structure around (7r/2,7r/2) rather following the nonin- 
teracting Fermi surface. It is especially hard to define a 
Fermi surface at all or extract a reliable estimate of the 
Fermi surface volume from these results. Note also the 
strong bifurcation around the X points. We find that es- 
pecially for small doping |Vn(k)| generally predicts that 
the Fermi surface bifurcates near X, with parts of the 
Fermi surface along both AI ^ X and X ^ T. This is 
in striking contrast to the fact that the spectra in Fig. |l| 
do not show a crossing along the direction X — > F, only 
a peak coming close to the Fermi energy. Apparently, 
|V?i(k)| is unable (at least at these temperatures) to dis- 
tinguish between a peak in the spectrum crossing or just 
simply approaching the Fermi energy. 

In accordance with the spectra in Fig. |[ the plot of 
|Vn(A:)| for large doping in Fig. |^ shows a fairly well de- 
fined Fermi surface that coincides with the Fermi surface 
of the non-interacting system. Again, in contrast to the 
sharp peaks found in the spectra, |Vn(fc)| shows a sub- 
stantial broadening, which in this case can however be 
explained by the conventional temperature broadening 
of Fermi's function. 



FIG. 4: |Vn(k)I for U = 2eV, T = -^eV, & = 0.20, = 
and Nc = 16. The color scheme is the same as in Fig. ^ 
The non-interacting and interacting Fermi surfaces coincide, 
in accord with Luttinger's theorem. 



An alternative way of mapping the Fermi surface, 
which also allows more direct contact with ARPES ex- 
periments, is to make constant energy plots of the single- 
particle spectra at the Fermi energy A(k, a; = 0). The 
constant energy plots A(k, — 0) are shown in Figs. ||, 
^ and for n = 0.95, n = 0.90 and n — 0.80, respec- 
tively, for the upper right quadrant of the first Brillouin 
zone. The regions of high density are colored in violet and 
low density in red. The solid black lines as before rep- 
resent the non-interacting Fermi surface. For n = 0.95, 
the Fermi surface resulting from our calculations is hole- 
like, centered around (tt, tt), and encloses a volume larger 
than the non-interacting Fermi surface, indicating a vi- 
olation of Luttinger's theorem. In addition, the shape 
of the Fermi surface close to (7r/2,7r/2), especially the 
clear shift above (7r/2,7r/2) cannot be interpreted nei- 
ther in terms of a simple tight-binding band structure 
nor a weak-coupling theorynS. Note also that in com- 
parison to the |V7i(k)| result, no apparent electron like 
shadow band at (tt, 0) can be seen. 

With increasing doping a bifurcation of the Fermi sur- 
face around (tt, 0) starts to develop (see Fig. signalling 
the incipient crossover from a hole like Fermi surface 
to the expected electron like at large doping. In addi- 
tion, around (7r/2,7r/2) the peak in the spectrum follows 
more or less the noninteracting Fermi surface again. This 
points towards a restoration of Luttinger's theorem. For 
n = 0.80, the Fermi surface is definitely electron like, cen- 
tered at (0,0), and has essentially the same volume and 
shape as the non-interacting surface, indicating a return 



FIG. 5: Constant energy scans of A{]i.,uj — 0) for U — 2eV, 
T = ^eV, S = 0.05, t± =0 and Nc = 16 in the right upper 
quadrant of the first BZ. The violet (red), represents regions 
of high (low) electronic density. The blue and violet regions 
map out the Fermi surface. The non-interacting Fermi surface 
is represented by a black line. The interacting Fermi surface 
is hole-like, centered at (7r,7r), and encloses significantly more 
volume than the non-interacting Fermi surface. This indicates 
a violation of Luttinger's theorem. 

to Fermi-liquid like behavior. No apparent remnants of 
the hole- like Fermi surface at lower doping are left. In- 
terpreting these bifurcations seen most strongly around 
n = 0.90 as shadow Fermi surfaces, we find a strong re- 
duction of the weight in these shadow features both for 
dopings less and larger than this "optimal" doping. In- 
terestingly, a similar behavior was recently observed in 
the analysis of ARPES data for Bi2212El. However, in 
these experiments the Fermi surface remained hole-like 
throughout the whole doping regime studied. Whether 
this might be related to a finite t' neglected in our present 
calculations will be discussed elsewhere. 



Summary and conclusions 

The increasing precision and quality of experimental 
ARPES spectra in recent years has led to a number of 
new results on the single-particle dynamics of the high- 
Tc cuprates, both partially resolving long-standing is- 
sues and posing new questions and problems. Motivated 
by especially the interesting observations concerning the 
changes of Fermi surface topology with doping, we have 
studied the two-dimensional (2D) Hubbard model with 
nearest neighbor hopping t in the intermediate coupling 
regime (on-site correlation U equal to the bandwidth) at 



FIG. 6: Constant energy scans of A(k,aj — 0) for 5 = 0.1, 
The other parameters and the color scheme are the same as 
in Fig. ^ The non-interacting and interacting Fermi surfaces 
start to coincide again, which can be interpreted as restoration 
of Luttinger's theorem. Note however the strong bifurcation 
around the X point. 



the temperature T — 0.033eV. To study the model in this 
most problematic parameter regime we used quantum 
Monte Carlo within the dynamical cluster approximation 
for the cluster size Nc = 16. Since this method allows 
for controlled and reliable calculations of low-energy fea- 
tures within a non-perturbative scheme and in the ther- 
modynamic limit, fundamental problems in this field can 
be addressed. These include the single-particle spectral 
properties at low energies, especially possible deviations 
from Luttinger's theorem or the formation of non Fermi 
liquid states, and the resulting topology of the Fermi sur- 
face 

^From the two different ways to define the Fermi sur- 
face, i.e. via |Vn(fc)| and inspection of a constant energy 
scan A{k,uj — 0) the latter turned out to be the more 
precise. The constant energy scans are able to distin- 
guish the situation where a peak in the spectral crosses 
the Fermi surface from that where it only approaches 
it. Thus they were free of spurious bifurcations observed 
in the |Vri,(fc)| plots. From the constant energy scans 
of the spectrum at the Fermi energy we find that the 
Fermi surface changes its topology compared to the non- 
interacting one as a function of decreasing doping. While 
the latter is closed around the zone center for every fi- 
nite doping, the interacting Fermi surface at low doping, 
6 — 0.05, is hole-like, closed around M = (tt, tt). More- 
over, the form and shifts present in the Fermi surface 
must be taken as clear evidence for a violation of Lut- 



FIG. 7: Constant energy scans of A(k,aj — 0) for 5 = 0.2. 
The other parameters and the color scheme are the same as 
in Fig. 1^ The non-interacting and interacting Fermi surfaces 
coincide, in accord with Luttinger's theorem. 



finger's theorem. Non Fermi liquid behavior is also evi- 
denced by the spectrum. A strong particle-hole asymme- 
try is found at the Fermi surface crossing near {tt/2, 7r/2) 
with hole-like excitations having much longer life-times 
than electronic excitations. Furthermore the correspond- 
ing self-energy shows strong /c-dependence rendering lo- 
cal approximations like the DMFA irrelevant. Additional 
structures in the self energy and a rather large residual 
scattering rate can also be interpreted as signs for non 
Fermi liquid behavior. However, to really distinguish a 
non Fermi liquid from a Fermi liquid with a possibly ex- 
tremely small energy scale much lower temperatures must 
be studied. 

With increasing doping the Fermi surface bifurcates 
at (tt, 0) around a doping 6 — 0.10 and from the topol- 
ogy around (7r/2,7r/2) one can infer a tendency to- 
wards restoration of Luttinger's theorem and conven- 
tional Fermi liquid behavior. Finally when 6 = 0.20 the 
Fermi surface is electron-like again, i.e. closed around 
the zone center, and coincides with the non-interacting 
Fermi surface. The corresponding spectra show well de- 
fined and sharp quasi-particle peaks around the Fermi 
energy with particle-hole symmetry being recovered in 
the low-energy region. 

Although the Hubbard model surely presents an over- 
simplification of the real cuprates, the common belief is 
that at least the essential qualitative features of the low 
energy dynamics should be reproduced. In particular, 
the results presented here show some interesting qualita- 
tive agreement with experiments regarding the behavior 
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of the different structures observed in the spectra. More- 
over, these features could be related to at least an ap- 
parent violation of Luttinger's theorem or possibly even 
non Fermi liquid behavior at low doping. From our re- 
sults it becomes very clear that neither weak-coupling 
treatments nor local theories like the DMFT are able to 
capture the essentials of the physics of the cuprate in 
the weakly doped regime. There are, of course, a vari- 
ety of further questions to be addressed. For example, 
most of the interesting and controversial ARPES results 
are for the system Bi2212, which actually is a bilayer sys- 
tem. Moreover, the simple nearest-neighbor tight binding 
band structure used in this paper is definitely not suffi- 
cient to describe the cuprates. Thus, a detailed study 
about the influence of a finite and negative t' and also 
a bilayer coupling on the topology of the Fermi surface 
or more generally on the low-energy single particle dy- 
namics is definitely necessary. We believe that such an 
investigation can also address some of the still mysteri- 
ous features in behavior of the Fermi surface topology 
of the high-Tc cuprates. Moreover it is very important 



to find new methods to solve the DCA self consistency 
cycle at lower temperatures or preferably at T = 0. This 
would enable us to clearly distinguish between a strong 
coupling Fermi liquid and genuine non Fermi liquid be- 
havior at low doping. 
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